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Abstract 

Diffusive transport of particles or, more generally, 
small objects is a ubiquitous feature of physical and 
chemical reaction systems. In configurations contain- 
ing confining walls or constrictions transport is con- 
trolled both by the fluctuation statistics of the jitter- 
ing objects and the phase space available to their dy- 
namics. Consequently, the study of transport at the 
macro- and nanoscales must address both Brownian 
motion and entropic effects. With this survey we re- 
port on recent advances in the theoretical and numer- 
ical investigation of stochastic transport occurring ei- 
ther in micro-sized geometries of varying cross section 
or in narrow channels wherein the diffusing particles 
are hindered from passing each other (single file dif- 
fusion). For particles undergoing biased diffusion in 
static suspension media enclosed by confining geome- 
tries, transport exhibits intriguing features such as (i) 
a decrease of nonlinear mobility with increasing tem- 
perature or, also, (ii) a broad excess peak of the ef- 
fective diffusion above the free diffusion limit. These 
paradoxical aspects can be understood in terms of 
entropic contributions resulting from the restricted 
dynamics in phase space. If, in addition, the suspen- 
sion medium is subjected to external, time-dependent 
forcing, rectification or segregation of the diffusing 
Brownian particles becomes possible. Likewise, the 



diffusion in very narrow, spatially modulated chan- 
nels gets modified via contact particle-particle inter- 
actions, which induce anomalous sub-diffusion. The 
effective sub-diffusion constant for a driven single file 
also develops a resonance-like structure as a function 
of the confining coupling constant. 

Keywords: diffusive transport, Brownian motion, en- 
tropic effects, diffusion enhancement, single file diffu- 



1 Introduction 

Effective control of transport in artificial micro- and 
nano-structures requires a deep understanding of the 
diffusive mechanisms involving small objects and, in 
this regard, an operative measure to gauge the role 
of fluctuations. Typical such situations are encoun- 
tered when studying the transport of particles in bi- 
ological cells [T] and in zeolites [SJ, catalytic reac- 
tions occurring on templates or in porous media [3], 
chromatography or, more generally, separation tech- 
niques of size disperse particles on micro- or even 
nanoscales [4]. To study these transport phenomena 
is in many respects equivalent to studying geometri- 
cally constrained Brownian dynamics [6l [7] . The 
fact that the diffusion equation is closely related to 
the time evolution of the probability density P(x,t) 
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to find a jittering particle at a location x at time 
t dates back to Einstein's pioneering work on the 
molecular-kinetic description of suspended particles 

With this minireview we focus on the problem 
of the diffusion of small size particles in confined 
geometries. Restricting the volume of the phase 
space available to the diffusing particles by means 
of confining boundaries or obstacles originates re- 
markable entropic effects .9,. The driven transport 
of charged particles across bottlenecks, such as ion 
transport through artificial nanopores or artificial 
ion pumps [TUl [Til H2 H3] or in biological chan- 
nels [H [HI OH UH US] is the better known model 
where diffusion is determined by entropic barriers. 
Similarly, the operation of artificial Brownian mo- 
tors [191 I5UI I5T] , molecular motors [55] and molecular 
machines [23] also results from the interplay of dif- 
fusion and binding action by energetic or, more rel- 
evant in the present context, entropic barriers [23] . 
The efficiency of such nanodevices crucially depends 
on the fluctuation characteristics of the relevant de- 
grees of freedom 25 . In addition, the interplay of 
diffusion over entropic barriers and unbiased time 
periodic drives is responsible for certain paradoxi- 
cal transport effects, like the recent observation of 
entropic-diffusion controlled absolute negative mobil- 
ity [2J. 

Although we restricted ourselves to small size par- 
ticle, we remind the reader that entropic forces surely 
affect the dynamics of extended chains diffusing on 
a periodic two (2D) or three dimensional (3D) sub- 
strate. A well established example is represented by 
the phonon damping of propagating solitons [27] . An- 
other example of great interest for its potential ap- 
plications in nanotechnology is the translocation of a 
long polymer molecule through a pore with opening 
size comparable with the polymer gyration radius. In 
this case, entropic effects were first predicted to theo- 
retically describe the conformation changes the chain 
undergoes to move past a conduit constriction [55J 
and, then, experimentally observed both in biologi- 
cal [55] and artificial channels [3D]. 

Another instance of constrained Brownian dynam- 
ics that rests, indeed, within the scopes of our minire- 
view, is the so-called single file diffusion. The motion 



of an assembly of small size particles in a narrow 
channel can be so tightly restricted in the transverse 
directions that the particles arrange themselves into 
a single file. The longitudinal motion of each par- 
ticle is thus hindered by the presence of its neigh- 
bors, which act as nonpassing movable obstacles. As 
a consequence, interparticle interactions in one di- 
mension can suppress Brownian diffusion and lead to 
the emergence of a new subdiffusive dynamics |31j . 

The outline of this work is as follows. In Sec. [2] 
we detail the geometry of the channels relevant to 
our review and set up the mathematical formalism 
needed to model the diffusion of a Brownian particle 
immersed in a confined suspension fluid. In Sec. [3] we 
introduce some exact results for the mobility and the 
diffusion coefficient of a driven Brownian particle in a 
one dimensional (ID) periodic potential, results that 
will be handy in the subsequent sections. In Sec. 2] 
we compute explicitly the entropic effects on particle 
transport in a static fluid filling a periodically mod- 
ulated channel. In Sec. [5] we address particle trans- 
port in a suspension fluid flowing along the channel 
subject to stationary pumping. In Sec. [6] we review 
numerical and analytical results for the diffusion of 
a single file along a periodically corrugated channel 
both in the presence and in the absence of an external 
drive. 

2 Channel models 

We consider the diffusive dynamics of spherical par- 
ticles in 2D or 3D pores, or channels, which extend 
in the x direction and possess a periodically vary- 
ing cross section. These channels are supposed to be 
symmetric with respect to a reflection on the channel 
axis in 2D (like those sketched in Fig. [TJ, or to any 
rotation about the channel axis in 3D. The channels 
are assumed to be delimited by rigid, smooth walls. 
The half width of the channel is described by the well- 
behaved boundary function w(x). The channel walls 
do confine the particles inside the channel but, oth- 
erwise, do not exchange energy with them. In partic- 
ular, we do not consider the possibility of adsorption 
of particles at the walls. Since we will disregard any 
rotatory motion of the particles about their centers of 
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Figure 1: (Color online) Brownian particles in a nar- 
row cylindrical channel directed along the x axis and 
with periodically modulated boundary w(x) (longitu- 
dinal section), (a) Pointlike particle suspended in a 
static fluid and subjected to a constant driving force 
F (Sec. 0|; (b) Single spherical particle dragged along 
by a laminar flow (Sec. O; (c) Single file diffusion of 
nonpassing driven particles (Sec. 16 . 2[) . 

mass, we need not specify particle- wall contact forces. 
Moreover, the radius of the particles is supposed to 
be smaller than the minima of w(x) (channel bottle- 
necks) , so that the particles are not restricted to stay 
within a confined region of the channel, but rather 
may diffuse everywhere along the channel. 

2.1 Diffusion equations 

The motion of particles that are immersed in a fluid 
medium may be influenced by various types of forces. 
We will here restrict ourselves to the discussion of 
small radius (almost pointlike) particles, whose pres- 
ence does not significantly modify the free motion of 
the fluid around them (laminar flow regime). The 
systematic impact of the fluid on the motion of par- 
ticles at position x = (x,y, z) is given by the Stokes 
force [32] 

^Stokes = t - V(x)] (1) 

where v(x, t) is the instantaneous velocity of the fluid 
in absence of the particle, x is the instantaneous par- 



ticle velocity, and 

7 = 6nrjR (2) 

is the friction constant, which is determined by the 
shear viscosity rj of the fluid, and the radius R of 
the particle. At the same time the fluid exerts on 
the particle also a random thermal force Fth- In the 
following, we only deal the case of fluids with homo- 
geneous temperature T and velocity fields that are 
almost constant on the particle scale R. As a conse- 
quence, to insure thermalization at temperature T it 
suffices to set 

F th (t) = (3) 

where fee denotes the Boltzmann constant and £(t) 
is a standard 3D Gaussian noise with (£(t)) = and 
(&(*)&(*')) = Sijd(t-t') for i,j = x,y,z. Other 
forces acting on the particles, like hydrodynamic in- 
teractions among different particles and between sin- 
gle particles and the wall, will be neglected [33J GH] ■ 
This simplification, in particular, requires a suffi- 
ciently low particle density. 

Finally, an external force F ext may act on the par- 
ticles describing e.g. the gravity force, or in the case 
of charged particles, an electrostatic force. We fur- 
ther specialize our discussion to constant longitudinal 
external forces pointing in the direction of the sym- 
metry axis of the channel. The dynamics of the center 
of mass x(t) of a single particle is then governed by 
Newton's equation of motion 

mi = F cxt - 7 [f - v(x, t)] + y/TrfkaT^t), (4) 

where we explicitly allow for a possible time- 
dependence of the fluid velocity field. For micropar- 
ticles moving with typical velocities of the order of 
1 cm/s, the inertial term mx(t) in Eq. ([4]) is neg- 
ligibly small compared to the environmental forces 
[35] : therefore, provided that the fluid velocity does 
not changes too fast, that is for spectral frequencies 
less than a few 100 Hz, one can safely set m = 
(overdamped limit or Smoluchowski approximation). 
Under these conditions the equation of motion (|4]) 
can be simplified as 



x = v(x,t) + ~F ext + J—?-Z(t). (5) 
7 VI 
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This Langevin equation is equivalent to the following 
Fokker-Planck equation [36] for the probability den- 
sity P(x, t) of a particle to be found at the position 
x at time t: 

^ = -V- fat). (6) 

Here J(x, t) denotes the corresponding probability 
current density 

J(x, t) = -( v(x, t) + — ) P(x, t) + ^VP(£, i). 
\ 7 J 7 

(7) 

These equations have to be supplemented by appro- 
priate boundary conditions, which will be discussed 
in the following subsection. 

2.2 Boundary conditions 

The channel is typically characterized by two bound- 
ary regions: A transverse boundary naturally results 
from the presence of the walls, whereas a longitudi- 
nal boundary is required to account for the channel 
length. 

The probability flux normal to the boundary in the 
presence of a rigid hard wall must vanish. Thus, to 
prevent pointlike particles from leaving the channel 
or being adsorbed at the walls, we must impose that 

n(x) ■ J{x, t) = x € wall, (8) 

where n(x) denotes the unit vector normal to the wall 
at point x. Note, however, that the center of mass of 
a spherical particle of finite radius R may approach 
the wall only up to a distance R, so that for finite 
size particles the boundary conditions (jSJ) applies on 
an appropriate inner surface parallel to the channel 
walls [37] . 

Various kinds of boundary conditions exist that 
regulate the inward and outward probability flows 
at the ends of a channel [32] . If the channel con- 
nects large, well mixed particle reservoirs, then con- 
stant probability densities Pl,r may be assigned at 
the channel ends. This leads to Dirichlet boundary 
conditions of the form 

P(x L ,y,z) = P Ll P(x R ,y,z) = P R , (9) 



where xl and x R denote respectively the left and 
right endpoints of the channel. A more detailed de- 
scription of the particle flow into and out of a channel 
can be achieved by relating flux and probability den- 
sities at x — xr.l, that is [35] 

J x {x L ,y,z,t) = -n L P(x Ll y,z,t) , 
J x (x R ,y,z,t) = K R P(x R ,y,z,t) . (10) 

Here, positive (negative) constants kl R correspond 
to partially absorbing (emitting) boundaries. As spe- 
cial cases, reflecting and absorbing boundaries corre- 
spond to k R l = and k Ri l = oo, respectively. 

Finally, for an infinitely long channel the periodic 
boundary conditions 

P{x,y,z,t)=P(x + L,y,z,t) (11) 

are more appropriate [36j . In the case of velocity 
fields which are constant with respect to time or vary 
periodically in time, these boundary conditions allow 
for stationary flux carrying solutions [T9j [20l [H] . 

3 Exact results for ID systems 

In order to set the stage, we consider first the ideal 
case where the diffusion of a particle in a periodi- 
cally corrugated channel can be assimilated to the 
diffusion on an energetic landscape represented, for 
simplicity, by a ID periodic substrate V(x) with pe- 
riod L, namely V(x + L) = V(x). Such a model is 
often employed to model, for instance, nanotube [39] 
and zeolite diffusion [2J. Let us consider a Brownian 
particle with mass m, coordinate x, and friction co- 
efficient 7, subjected to a static external force F and 
a thermal noise F t h(t). In the notation of Sec. 12.1 [ 
we set F cxt = -V'{x) + F. 

The corresponding stochastic dynamics is de- 
scribed by the Langevin equation 

mi = -V'(x) -jx + F+ v/2^fef £(t), (12) 

where £ (t) is the standard Gaussian noise also defined 
in Sec. 12.11 Moreover, for the purposes of this review, 
the substrate V(x) can be taken symmetric under 
reflection, V(x) = V(—x). 
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In extremely small systems, particle fluctuations 
are often described to a good approximation by the 
overdamped limit of Eq. (|12[) . i.e., by the massless 
Langevin equation 



and for the diffusion coefRcient 



7± = -V'(x) +F+ , 



(13) 



where the inertia term mx has been dropped alto- 
gether (Smoluchowski approximation). 

An overdamped particle is trapped most of the 
time at a local minimum of the tilted substrate as 
long as F < F3, F3 denoting the depinning threshold 
F3 = m&x{V'(x)}. Drift occurs by rare noise in- 
duced hopping events between adjacent minima. For 
F > F3 there exist no such minima and the parti- 
cle runs in the F direction with average speed ap- 
proaching the no-substrate limit F/j. This behavior 
is described quantitatively by the mobility formula: 



F 



with 



lim 

t — >oc 



(*(*)) 

t 



L 



(14) 



(15) 



(t(L,F))- 

Here and in the following, (t n (L,F)} denotes the n- 
th moment of the first passage time of the particle 
across a substrate unit cell in the F direction. 

As the particle drifts subjected to the external force 
F, the random hops cause a spatial dispersion of the 
particle around its average position (x(t)). The cor- 
responding normal diffusion coefficient, 



D(F) = lim 

t—>oo 



2t 



(16) 



can be computed analytically by regarding the hop- 
ping events in the overdamped regime as manifesta- 
tions of a renewal process, that is [40] : 



D(F) = L 



2 {t 2 {L,F))-{t{L,F)Y 



2{t(L,F)Y 



(17) 



Simple algebraic manipulations lead to explicit ex- 
pressions for the nonlinear mobility 36 , 4TJ [42l [43] , 



DqL 1 



F 



-LF/k B T 



Jo I+{x)dx 



(18) 



D(Fl = L 2 Jo I 2 +{x)I^{x)dx 
Do [J L I + {x)dxY 



(19) 



Here, I ± (x) = ^ e^ ±v WTV{xTy)-yF}/k B T dy and 
^0 = k-gT/j denotes Einstein's coefficient for a free 
diffusing Brownian particle. 

For F — > Eqs. (JTHJ) and (TJ5]) reproduce the 
zero-bias identity D(0)/D o = 7A«(0) with A*(0) = 
L 2 / (jJ L I+(x)dx)\M H3- Notably, as F ap- 
proaches the depinning threshold F3 the mobility 
curve (|18p jumps from zero (locked state) up to close 
I/7 (running state). Correspondingly, the diffusion 
coefficient (IT9l) develops a diffusion excess peak, i.e. 
with D > Do, consistently with numerical observa- 
tions in Ref. [55]. Both the mobility step and the 
D peak get sharper and sharper as T is lowered [40) . 
The same conclusions apply in the presence of iner- 
tia (i.e. for particles of finite mass m) [45j [46], as 
well, with the difference that the relevant depinning 
threshold shifts towards lower values (proportionally 
to 7 as 7 — > [3"6]1. 

Finally, we emphasize that the above formulas 
for the nonlinear mobility and the effective diffu- 
sion coefficient retain their analytic structure also 
when generalized to anomalous (sub)-diffusion on a 
ID substrate, by merely substituting the normal dif- 
fusion constant Do by the fractional diffusion con- 
stant occurring in the corresponding fractional diffu- 
sion equation [4"T] . 



4 Particle diffusion in a static 
fluid 

In the limiting case of pointlike particles with zero 
interaction radius diffusing in constrained 2D or 3D 
geometries, elastic contact particle-particle interac- 
tions can be neglected. As long as hydrodynamically 
mediated interactions can also be neglected, see e.g. 
[33l [34] , the confining action of the channel walls is 
modeled by the perfectly reflecting boundary condi- 
tions (El). 
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In the presence of a constant external force point- 
ing in the channel direction, an overdamped Brown- 
ian particle suspended in a static medium is described 
by the Langevin equation (J5]), or by the correspond- 
ing Fokker-Planck equation ©, both with zero ve- 
locity field v(x, t) = 0. For a general choice of the 
periodic boundary w(x), there exist no exact analyt- 
ical solutions to the Fokker-Planck equation © with 
boundary conditions l[5]). However, approximate so- 
lutions can be obtained by reducing the problem of 
free Brownian diffusion in a 2D or 3D channel to that 
of Brownian diffusion on an effective ID periodic sub- 
strate (Sec. [3]). Under such a scheme, narrow channel 
constrictions, corresponding to geometric hindrances 
in the fully dimensional problem, are modeled as en- 
tropic ID barriers. 

In the absence of external forces, i.e., for F ext — 
0, particle dynamics in confined structures [see 
Fig. [Ha)] can be approximately described by the 
Fick-Jacobs kinetic equation with spatially depen- 
dent diffusion coefficient gi M GHH E2 : 



dP{x,t) 
dt 



d_ 

dx 



D(x)a(x) 



d P 



dx a(x) 



(20) 



with a{x) denoting the dimcnsionless channel cross- 
section 2w(x)/L in 2D and irw 2 (x)/L 2 in 3D. Equa- 
tion PD|) was obtained [35] from the full Fokker- 
Planck equation ^ for small amplitude boundary 
modulations w(x), on assuming a uniform density 
P(x, t) and integrating out the transverse coordinates 
[e.g. for a 2D channel P(x,t) = dy P(x, y, t)]. 

At variance with the original Fick-Jacobs equation 
[48j . introducing an x-dependent diffusion coefficient 
considerably improved the accuracy of the kinetic 
equation (|20p thus extending its validity to larger 
w(x) amplitudes [49], [50l [51] . Here, the expression 



D(x) 



Dn 



[1 + w'(x) 



2]a 



(21) 



with a = 1/3 in 2D and a = 1/2 in 3D, was de- 
termined to best account for wall curvature effects 

In the presence of weak longitudinal drives F, the 
Fick-Jacobs equation ([2H)l can be further extended to 



[50115211531: 



™ = ±D(x) (™ + ^P 
8t dx { '[dx + k B T 



(22) 



where the free energy A(x) = E{x) — TS(x) is made 
up of an energy, E{x) = —Fx, and an entropic term, 
S(x) = k B \o.a{x). For a periodic channel, A(x) as- 
sumes the form of a tilted periodic potential. More- 
over, for a straight channel, w'(x) = 0, the entropic 
contribution vanishes altogether and the particle is 
subject to the sole external drive. Of course, for 
F = the free energy is purely entropic and Eq. (|22| 
reduces to the Fick-Jacobs equation (|20p. An alterna- 
tive reduction scheme based on macrotransport the- 
ory has been proposed recently in Ref. [5^] . 

Key quantifiers of the reduced ID kinetics (|2"2"|) arc 
the average particle current or, equivalently, the non- 
linear mobility ([T4]) . and the effective diffusion coef- 
ficient (fT7|) defined in Sec. [U On generalizing the 
derivation of Eqs. (fT8|) and (fT9|) to account for the 
x-dependence of D(x), one obtains [53], respectively, 



7M(/) 



1 



-/ 



(23) 



and 



Pes 
Do 



dx 



dz D(z) e 



A(x)/k B T 



[i{z)Y 



where 



I(x) 



a A(x)/k B T px 



D(x)/D 



dy _ 



A(y)/k B T 



(24) 



(25) 



and the dimensionless force / is defined as the ratio 
of the work done by F to drag the particle a distance 
L, to the thermal energy fceT 1 , that is 

FL 



f 



k B T 



(26) 



We stress here an important difference with the 
energetic ID model of Sec. [3] For a particle mov- 
ing in a ID periodic substrate V(x), the barriers AV 
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separating the potential wells provide an additional 
energy scale besides FL and k#T, so that the particle 
dynamics is governed by at least two dimensionless 
energy parameters, say, AV/ksT and FLjk^T [55] . 
In contrast, Brownian transport in a periodically cor- 
rugated 2D or 3D channel is solely determined by 
the dimensionless force / [53] . This can be proven 
by rescaling the problem units as follows. We mea- 
sure all lengths in units of the period L and time in 
units of t = L 2 / Dq, that is, twice the time the parti- 
cle takes to diffuse a distance L at temperature T. In 
these new dimensionless variables the Langevin equa- 
tion ([5]) for a Brownian particle in a static medium 
reduces to the equation 



*(*) = / + m , 



(27) 



with no tunable constants [55] , As a consequence, 
the ensuing Brownian dynamics is controlled only by 
the parameter / and, of course, by the additional 
length(s) that possibly enter the boundary function 
w(x). 

The /-dependence of the average particle current 
(|23| and the effective diffusion ([24]) were compared 
with the results obtained by numerical integration of 
the 2D Langevin equation J5]) (for details see Ref. 
[52|). For simplicity, the channel walls were assumed 
to have the sinusoidal profile 



w(x) = a[sin(27ra;/L) + k] 



(28) 



with a > and k > 1. Here, a(n ± 1) are, re- 
spectively, the maximum and the minimum chan- 
nel width. Moreover, a/L controls the slope of the 
boundary function w(x), which in turn determines 
the modulation amplitude of the diffusion coefficient 
D(x) in Eq. plj). The mobility fi, Eq. (fH|). and 
the corresponding effective diffusion coefficient D, 
Eq. (|f 6p . were computed by ensemble averaging 3T0 4 
simulated trajectories. 

The most significant results are displayed in Figs. [5] 
and[3J At variance with the purely energetic ID mod- 
els of Sec. [31 the nonlinear mobility decreases upon 
increasing the strength of thermal noise. Moreover, 
an enhancement of the effective diffusion coefficient, 
with maximum exceeding the free diffusion constant 
Dq, was also observed [53] 155] . 
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Figure 2: (Color online) Nonlinear mobility /i vs. 
driving force in dimensionless units, / = FL/kBT, 
for a 2D channel at different temperatures k^T = 
0.01 (crosses), 0.1 (pluses), 0.2 (squares), and 0.4 (tri- 
angles). After rescaling, all data sets collapse on one 
curve which at low / closely compares with the an- 
alytic approximation (|23|) (solid curve). Other sim- 
ulation parameters are L = 1, 7 = 1, and w(x) — 
[sm(27nr) + 1.02] /2tt. The inset shows the devia- 
tion of the analytic approximation (|23| (solid curve) 
from the numerical results (dotted curve) for large /. 
Note that the numerical curve approaches the correct 
asymptotic limit 7/x = 1. 



At low values of the control parameter / the ana- 
lytical approximations ([23)) and ([24)) match perfectly 
the corresponding numerical curves, whereas devia- 
tions occur at high /. Most remarkably, contrary 
to the simulation output, the analytical curves for 
D e g I Dq and 7/x fail to attain the correct asymptotic 
limit 1 for / — * 00 (see insets in Figs. [5] and [3]). This 
occurs because the assumption of uniform density dis- 
tribution, introduced in the Fick- Jacobs formalism to 
eliminate the transverse coordinates, is no more ten- 
able in the presence of strong drives. 

The agreement between theory and numerics im- 
proves for smooth modulations of the channel walls, 
i.e. for small boundary slopes |u/(:r)|, which can be 
achieved, e.g., for appropriately small a [5U[SS]. A 
phenomenological criterion to asses the validity of the 
stationary state solutions of the Fick- Jacobs equation 
([22]) can be formulated by comparing the different 
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Figure 3: (Color online) Effective diffusion coefficient 
-D e ff vs. / for the same simulation parameters as in 
Fig. [2 Here, too, the rescaled data collapse on one 
curve, which asymptotically approaches the correct 
limit D e g / Dq = 1 (inset). The analytic approxima- 
tion ([24]) (solid curve) fits well the raising branch of 
the numerical data sets. 

characteristic time scales of the problem, namely 
(i) the diffusion times, respectively, in the transverse, 



(2a 2 /A))(l + K) 2 



(29) 



(ii) and in the longitudinal direction, 

rj d) = L 2 /2D a , (30) 

(iii) and the drift time the applied force F takes to 
drag the particle one channel unit length L across, 



rP= 7 L/F. 



(31) 



Uniform probability distribution in the transverse 
direction sets on only if the transverse diffusion mo- 
tion of the particle is sufficiently fast relative to 
both the diffusive and the drift longitudinal motions, 
which implies t± <C min{r,j , r,| ■*}. For large drives 

it suffices to require that -C T \\ \ which leads to 
the condition 1551: 



f<fc 



1 



2(1 + Kf 



L 



(32) 



100 



50 




Figure 4: Critical force f c vs. boundary modulation 
amplitude a for the 2D channel of Fig. [2j Other sim- 
ulation parameters are L = 1, 7 = 1, and k^T = 0.1. 
The numerical values of f c (solid dots) were obtained 
by requiring deviations of Fick- Jacobs approximation 
from the simulation results smaller than ~ 1%. The 
solid line is the fitting law f c = Ci(L/a) 2 + C2 with 
c\ = 0.792 and C2 = 3.524. The apparent discrep- 
ancy with Eq. (f3"2"]> results from the above-mentioned 
numerical accuracy criterion and is expected to van- 
ish for higher accuracy. Inset: same data sets on a 
logarithmic scale. 



The critical force parameter f c , above which the 
Fick-Jacobs description is expected to fail, depends 
on the remaining free parameters of the problem. For 
the boundary function (|28p . fc is a function of a/L, 
alone, as illustrated in Fig. [4] For a L the Fick- 
Jacobs approximation becomes untenable already for 
relatively small forces /, whereas for a <C L its valid- 
ity extends to significantly larger drives. 

5 Particle diffusion in moving 
fluids 

Let us assume now that a small, spherical particle is 
swept along in the stationary velocity field w(x, t) of 
a moving, incompressible fluid [see Fig.[ljb)], rather 
than by a constant force, F cxt , like in Sec. [4j For 
-Fext = 0, its dynamics is then described by the 
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Particle diameter [p.m] 

Figure 5: Average induced particle current v G ver- 
sus particle diameter. The particle is immersed in 
a laminar flow confined to a channel with broken 
reflection symmetry (inset). The fluid, with viscos- 
ity vr = rj/Vwatex and 77 watC r = 1-025 • 10 3 Ns/m , is 
pumped periodically back and forth, along the x axis 
of the channel, with frequency v. For further details, 
see Ref. [56] . 

Langevin equation, in the overdamped limit, 

S(t) = v(x(t), t) + sj^^m , (33) 

where the fluid velocity field v(x, t) is to be deter- 
mined. 

To this purpose, we write first the Stokes equa- 
tions for a stationary incompressible flow in the limit 
of vanishing Reynolds numbers (the convective accel- 
eration terms in the Navier-Stokes equations can be 
safely neglected). One obtains the so-called "creeping 
flow equation" [32] : r]Av(x) = Vp(x) , where p(x) is 
the pressure field responsible for the stationary lam- 
inar flow of the fluid and r\ is the absolute dynamic 
fluid (or shear) viscosity. On introducing the scalar 
field &(x,r), with (r, <f>) being the polar coordinates 
in the transverse plane (y,z), v(x) can be rewritten 
as 32 

v(x, r) = V x r)e /r) , (34) 
where e$ is the unit <f> vector. On substituting 



Eq. ([34")) into the creeping flow equation, one obtains 
the following linear homogeneous fourth-order differ- 
ential equation [3H [55] : 

(did d 2 \ 2 

( r *;* + s?) *m = °. (35) 

with boundary conditions 

*|r=0 = C, (36) 

1:^=0 = 1^*1-0 = 0, (37) 

VW(x,r = w(x)) =0, (38) 
V(x + L,r) = %(x,r), (39) 

where c is an arbitrary constant. Solutions to 
Eq. (|35[) . determine the velocity field v(x) up to a 
multiplicative factor, which, in turn, can be estab- 
lished by imposing some additional condition, e.g. for 
the pressure drop across a channel unit, i.e. [56j 

p(x,r)-p(x + L,r) = -2jf dx. (40) 

Once the velocity field for a particular channel ge- 
ometry and pressure profile is known, the Langevin 
equation (|3"3"j) can be solved numerically. 

Particle separation across a micro-channel is a pro- 
cess of great technological importance [7j [19]. Any in- 
homogeneity in the spatial distribution of an ensem- 
ble of non-interacting suspended particles can only 
be caused by the hydrodynamic interaction between 
particles and walls [37]. The no-flux boundary con- 
dition ([8]) for unbiased particles with F GXt = reads: 

n(x)- (v(x)P(x, t) - ^-VP(x, t)j = , X G wall . 

(41) 

In the case of vanishing drift velocity normal to the 
channel walls, n(x) ■ v{x) — 0, only uniform distri- 
butions are allowed and no particle separation could 
ever be achieved [37]. However, as anticipated in Sec. 
12.21 for finite size particles the no-flux boundary con- 
dition |4T|) strictly holds on an effective inner sur- 
face at a distance from the channel walls. Due to 
its finite size, a particle cannot move steadily along a 
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given flow streamline as this gets too close to the 
walls; upon hitting the wall, the particle bounces 
into inner flow streamlines as if they were subject 
to a transverse field gradient. These hydrodynamic 
forces may, indeed, lead to accumulation and deple- 
tion zones within the channel. 

By generalizing such a boundary effect, Kettner 
et al. [56] predicted by numerical simulation that 
a micron-sized channel with broken reflection sym- 
metry can be used to separate particles according to 
their size, as illustrated in Fig. \5\ A time oscillating 
pressure profile, p(x,t) = p(x)f(t), where f(t) is a 
sinusoidal function with frequency v, was assumed to 
control the periodic flow of the fluid, back and forth 
along an infinitely long channel. Within the creep- 
ing flow approximation, the ensuing time dependent 
velocity field v(x,t) — v(x)f(t) was then obtained in 
terms of the solution v(x) of the unperturbed Stokes 
equation introduced above, and the corresponding 
Langevin equation ()33j) numerically integrated. Later 
on, this mechanism was demonstrated experimentally 
by Matthias and Muller [57]. 



6 Single File Diffusion 

As the design and the operation of biology inspired 
nanodevices [19] have become experimentally afford- 
able, understanding particle diffusion in a ID sub- 
strate has been recognized as a key issue in trans- 
port control [35]. In this context the particle- and 
particle- wall interactions play a central role. Pair in- 
teraction between thermally diffusing particles does 
not affect the normal character of Brownian diffusion, 
as long as the particles are able to pass one another, 
no matter how closely they are confined. This holds 
true even when, under appropriate temperature and 
density conditions, attracting particles cluster or con- 
dense in the substrate wells [58] . 

Things change dramatically for strictly ID geome- 
tries. Let us consider, for instance, an ensemble of N 
unit-mass particles moving with preassigned dynam- 
ics along a segment of length /. If the interparticle 
interaction is hard-core (with zero radius), the elastic 
collisions between neighboring particles are nonpass- 
ing - meaning that the particles can be labeled ac- 
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Figure 6: (Color online) Velocity autocorrelation 
function C(t) in a sSF at k^T = 1 and for (p, 7) = 
(0.1,3) (top, orange curve), (0.1,5) (middle, black 
curve), and (0.05,5) (bottom, green curve). Inset: 
log-log plot of the negative C(t) tails; the solid line is 
the predicted tail -C{t)/{D p) 2 ~ (t / 't a )~ § / (40F) . 



cording to a fixed ordered sequence. The particles are 
thus arranged into a file where, at variance with the 
situation described in Sec. their diffusion is sup- 
pressed by the presence of two nonpassing neighbors, 
also movable in the longitudinal direction. 

The diffusion of a free single file (SF), i.e. in the 
absence of a substrate, has been investigated in de- 
tail 131, ■ In the thermodynamic limit (I, N — > 00 with 
constant density p = N/l) the mean square displace- 
ment of each file particle can be written as 



(Ax 2 (t)) = \Ax(t)\/p 



(42) 



with |Ax(i)| denoting the absolute mean displace- 
ment of a free particle. For a ballistic single file (bSF), 
clearly |Aa:(t)| = (\v\)t, where (...) is the ensemble 
average taken over the distribution of the initial ve- 
locities, and therefore 



(Ax 2 (t)) = (\v\)t/p. 



(43) 



A bSF particle diffuses apparently like a Brown- 
ian particle with normal diffusion coefficient Dq = 
(\v\)/(2p). For a stochastic single file (sSF) of Brow- 
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nian particles with damping constant 7 at tempera- 
ture T, the equality |Ax(i)| = y/AD^TpK yields the 
anomalous diffusion law 

(Ax 2 (t)) =2D SF Vt/p, (44) 

where the mobility factor Dsf = V Dq/tt is related 
to the single particle diffusion constant Dq — k^T/j. 
The onset of the subdiffusive regime ([44]) occurs for 
t > Td, Td — (Dop 2 ) -1 being the average time a single 
particle takes to diffuse against one of its neighbors 
[5"i*"l I6"0"] . The diffusive regimes (03} and (gU) have 
been demonstrated both numerically [61 and exper- 
imentally 34. 62J. Let x(t) represent the coordinate 
of one file particle, assumed to be a continuous differ- 
entiable stochastic process with (x(t)) — 0. Kubo's 
relation [53] 

1 -^{A 2 x{t)) = fc{j)dT (f-oo), (45) 

with C(t) = (x(t)x(0)), best illustrates the role of the 
dimensional constraint on SF diffusion. In the case of 
normal diffusion, (A 2 x(t)) = 2Dot, the r.h.s. of Eq. 
([4"5"]) converges to the positive value Dq = J Q C(r)dr, 
namely x{t) diffuses subject to Einstein's law. The 
subdiffusive dynamics ([44]) . instead, is characterized 
by J °° C(r)rfr = 0, which, in view of Eq. (j45|) . im- 
plies that C(t) develops a negative power-law tail 
C(t) cpt' 13 with (3 = § and cp = D SF /Ap. Nu- 
merical simulations support this conclusion. In Fig. [6] 
we report a few curves C(i) for different p and 7. The 
negative tails are apparent and compare well with the 
estimate derived from Eq. (|45l) (see inset). They are 
a typical signature of the collisional dynamics in a SF: 
As an effect of pair collisions, an initial velocity x(0) 
is likely to be compensated by a backflow velocity of 
opposite sign [501 IM] ■ 

6.1 Stochastic single file on a sub- 
strate 

We focus now on the case of a SF diffusing on a si- 
nusoidal substrate with potential [65j [66] 

V(x) = d[l - cos(2ttx/L)}. (46) 




d/k B T 

Figure 7: (Color online) Diffusion of a stochastic sin- 
gle file in the periodic potential ([4*6"]) [5B] : (a) (Ax 2 (t)} 
vs. t for k^T = 1, d = 1 and increasing 7 (from 
top to bottom on the l.h.s); (b) (Ax 2 (t)} vs. t for 
7 = 5 and increasing d (from top to bottom on the 
r.h.s). The t and t 1 ^ 2 slopes (dashed lines) have been 
drawn for reader's convenience; (c) the diffusion mo- 
bility Dsf vs. d/k-oT for 7 = 5 and d = 1 (circles) 
and k^T = 1 (squares). Other simulation parame- 
ters are: N = 3 x 10 3 , L = 2tt, / = 3 x 10 3 L, and all 
particles have unit mass. 

This variation of the SF model rises naturally in con- 
nection with quasi- ID situations where the particles 
(not necessarily suspended in a fluid) can be repre- 
sented by disks moving along a narrow spatially mod- 
ulated channel with cross-section smaller than twice 
the disk diameter [see Fig. [l][c)]. The confining ac- 
tion of the channel can be accounted for by modeling 
the SF dynamics in terms of a periodic substrate po- 
tential with a effective strength d. This is the case, 
for instance, of most nanotubes and zeolite pores [2]. 

In the simulations of Rcfs. [60, 66J the i-th particle 
was assigned: (i) random initial position, Xi(0), and 
velocity, Xj(0). Upon each elastic collision it switched 
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velocity with either neighbors without altering the file 
labeling; (ii) independent Brownian dynamics deter- 
mined by a viscous force —jii and a thermal force 
\/2'-fkBT£ i i(t), with £i(t) uncorrelated standard Gaus- 
sian noises defined as in Sec. 12. 1[ in order to guaran- 
tee thermalization at temperature T. 

Numerical evidence led to conclude that the pe- 
riodic substrate potential V(x) does not invalidate 
the sSF diffusion law (|44p. although the dependence 
of the mobility factor Dsf on the system parame- 
ters becomes more complicated. In panels (a) and 
(b) of Fig. [7] Z?sf is characterized as a function of 
7 and T, respectively. The identity Dsf = \J -Do/7, 
reported above for V(x) = 0, applies here, too, un- 
der the condition of replacing Dq with the modified 
diffusion constant D(F) of Eq. In Fig. E£a) the 

rescaled curves 7 1 / 2 (Ax 2 (t)) versus t overlap asymp- 
totically for any damping regime. 

The temperature dependence of the mobility 
Dsf = y/ D(F)/ , y, is more interesting. As implicit 
in the discussion of Sec. El this prediction gets more 
and more accurate for large 7 and increasingly high 
activation-to-thermal energy ratios d/ksT. As a con- 
sequence, the rescaled mobility y/j/d Dsf turns out 
to be a function of d/k^T, alone, in good agreement 
with the simulation results displayed in Fig. EJc). 

6.2 Driven single files 

Finally, we consider the case of a driven sSF, namely, 
we now assume that all file particles are subjected 
to an additional constant force F pointing, say, to 
the right (F > 0). We know from Sec. [3] that the 
diffusion of a single Brownian particle drifting down 
a tilted washboard potential exhibits enhanced nor- 
mal diffusion with diffusion constant (fT9|) . Extensive 
simulation of a driven SF [66j yielded the numerical 
data reported in Fig. [5] The kinetic mobility of the 
file, defined as /1 = (x) [F, turned out to coincide 
with the mobility of a single particle under the same 
dynamical conditions, irrespective of 7 (inset of Fig. 
[8|). More remarkably, the subdiffusive regime (|44l) 
applies also in the presence of bias, though with an 
D-dependent mobility factor. When plotted versus 
d, Dsf attained a maximum enhancement for d > F, 
i.e., in coincidence with the (noise-assisted) depinning 
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Figure 8: (Color online) Diffusion of a driven stochas- 
tic single file in the periodic potential flS)) : the mobil- 
ity factor £>sf vs d for k-gT = 0.3, 7 = 5, and F = 
(squares) and F = 1 (circles). The solid curves rep- 
resent the law D SF = ^D{F)/tt with D(F) given m 
Eq. (fl9l) . Other simulation details are as in Fig. [7] 
Inset: the kinetic mobility, /i = (x)/F vs. Fj\fd for 
7 = 0.5 (circles) and 7 = 2 (squares) at k^T = 0.1. 
The fitting curves are the analytical predictions for a 
single Brownian particle based, respectively, on Eqs. 
(11.194) (low 7) and (11.50) (high 7) of Ref. [36]. 

of the file from its sinusoidal substrate. Again, the 
identity D$p — \J F>{F) /it combined with formula 
(fl~9)) for Do provided an excellent fit of the simulation 
data for large 7. Of course, the mobility enhance- 
ment at depinning can be revealed also by plotting 
Dsf versus F at constant d. 

7 Conclusions and Outlook 

With this minireview we presented the state of the art 
of diffusive transport occurring in systems character- 
ized by confinement due to either their finite size or 
to particle interactions in restricted geometries. Con- 
finement plays a salient role in the Brownian motion 
of driven particles. Indeed, the entropic effects asso- 
ciated with confinement may give rise to anomalous 
transport features. Some main new phenomena are 
the entropy-driven decrease of mobility with increas- 
ing temperature, the diffusion excess above the free 
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diffusion limit and the resonant behavior of the ef- 
fective diffusion coefficient as a function of control 
parameters such as temperature, external gradients, 
or geometry design. 

The reviewed transport features can thus be im- 
plemented in the design of new transport control se- 
tups or protocols. In particular, Brownian transport 
through confined geometries is expected to find appli- 
cation in nanotechnology, for instance, to operate de- 
vices for the separation of nanoparticles, to speed up 
schemes for catalysis on templates, or also to realize 
efficient, driven through- flows of microsized agents, 
or reactants, in miniaturized lab-on-a-chip devices. 
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